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The exact nature of the QCD phase transition has still not been determined con- 
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scaling behavior for two quark flavors. Ultimately, this issue can be resolved only 
by a careful scaling and finite-size scaling analysis of the lattice results. We use a 
renormalization group approach to obtain finite-size scaling functions for the 0(4)- 
model, which are relevant for this analysis. Our results are applicable to Lattice 
QCD studies of the QCD phase boundary. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) at finite temperature and density is currently actively 
researched both on the experimental and the theoretical side. The equation of state of QCD 
and, in particular, the nature of the phase transition from the strongly interacting hadronic 
phase to the strongly interacting quark-gluon plasma phase is of great importance for a 
better understanding of the experimental data [1]. 

In QCD, two phase transitions take place at finite temperature and density: a deconfine- 
ment phase transition dominated by the gauge fields and a chiral phase transition driven by 
the interplay between quark and gauge-field degrees of freedom. The deconfinement transi- 
tion in pure SU(3) gauge theory is of first order, but becomes a crossover in the presence of 
dynamical quarks. The nature of the chiral phase transition, in turn, depends on the num- 
ber of quark flavors, the strength of the explicit chiral symmetry breaking, and the strength 
of the chiral anomaly [21 |3]. Whether the chiral and the deconfinement phase transition 
coincide is still under investigation. Assuming that the effects of the anomaly at the chiral 
phase boundary are small, one expects a second order phase transition for two quark flavors 
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in the chiral limit and that QCD falls into the 0(4) universality class |2]. In this case the 
phase transition in QCD would be dominated by the restoration of chiral symmetry. With 
explicit symmetry breaking, the order of the phase transition changes instantaneously from 
second order to a crossover. If the strength of the chiral anomaly is large, however, for two 
massless quark flavors one expects that the transition is of first order [21 [3] . 

For studying full QCD, lattice simulations are currently the most powerful tool. However, 
the determination of the nature of the chiral phase transition still remains a difficult task, 
since such simulations are necessarily performed in finite volumes and the implementation 
of chiral fermions is difficult. While there is much evidence that QCD with two flavors 
falls into the 0(4) universality class [4], results obtained with two dynamical flavors of 
staggered fermions do not exhibit the expected scaling behavior [5l El [7] or suggest a flrst- 
order transition [H |9] . 

The flnite simulation volume of lattice simulations poses in particular a problem for the 
investigation of phase transitions. Since phase transitions occur strictly speaking only in 
the inflnite-volume limit and a continuous symmetry cannot be spontaneously broken in a 
flnite volume, the introduction of explicit symmetry breaking in the form of a flnite current 
quark mass term is mandatory. This makes it difficult to determine the nature of the chiral 
phase transition from lattice QCD results. An important tool for the analysis of lattice QCD 
data is the investigation of the finite-size scaling behavior. The underlying universality class 
determines the scaling behavior of e.g. the order parameter characterizing the transition 
and the corresponding susceptibility. If universal behavior obtains, results are expected to 
fall onto universal scaling curves, characterized by critical exponents and scaling functions of 
the underlying universality class. Thus the order of the phase transition and its universality 
class can be established by a comparison with the known critical exponents and the known 
scaling behavior. From such an analysis, there is indeed evidence of 0(4) scaling for QCD 
with Wilson fermions [TT| \T2\ [T3] . however, the expected scaling behavior has not been 
seen with staggered fermions |5], El [7j. Results with a modified QCD action with two fiavors 
of staggered fermions suggest that current simulation volumes might be actually outside 
the finite-size scaling region [H]. This result underlines the importance of a finite- volume 
scaling analysis of 0{N) models. 

So far, scaling and finite-size scaling functions have been determined mainly from lattice 
simulations of 0{N) spin models [T51 [TBI El- These results have already been used in 
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the scaling analysis of Lattice QCD results [H [TTl [121 113 IH]- this paper, we provide 
the technical framework for finite-size scaling studies with the functional renormalization 
group (RG) and compute the finite-size scaling functions for the linear 0(4) model. Our 
RG approach complements the one taken in 0{N) lattice simulations. It is computationally 
efficient and allows us to study scaling over a wide range of volume sizes and values of the 
external symmetry breaking field, which is usually associated with an external magnetic 
field in spin models or, in the case of QCD, with the current quark mass. The results are 
directly applicable to a comparison with Lattice QCD data [TH]. Making use of universality 
arguments, our approach enables us to study the scaling of the order parameter and the 
associated susceptibilities at the phase boundary where the fermions are assumed to have 
decoupled from the critical fiuctuations. While the approach in the present paper cannot be 
used to study the onset of chiral symmetry breaking in terms of quark-gluon dynamics [20l 
2T| [22] . it still provides important information for an analysis of Lattice QCD data which 
helps to shed light on the ongoing discussion about the nature of the QCD phase transition. 
The paper is organized as follows: In Sect. |TT], we present the setup of our RG formalism. 



In Sects, in and IV, we discuss general aspects of finite- volume and finite-size scaling in 
quantum field theories. Our results are then presented in Sec. IVl Concluding remarks and 



future plans are given in Sect. VI 



II. RENORMALIZATION GROUP APPROACH FOR FINITE VOLUME 

STUDIES 

In this section we discuss our Renormalization Group approach to studying finite-volume 
scaling. In the first part of this section, we briefiy discuss the derivation of the fiow equations 
for the 0(4)-potential in infinite volume. A detailed discussion of the derivation and the 
underlying approximations can be found in Ref. |23]. In the second part of this section we 
generalize our fiow equations to finite- volume and recapitulate some earlier results [2ll [25], 
[26]. 
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A. 



Infinite volume 



The effective action of the 0(4) model in d = 3 spatial dimensions is given by 




(1) 



where 0"^ = (cr, vf). The potential U{a,7i^) depends on a and vf^ separately, since the 
presence of a no n- vanishing external source term {—Ha) in the ansatz for the effective 
action is indispensable for a study of phase transitions in finite volume. We study the 0(4) 
model in the so-called local potential approximation (LPA), where we neglect a possible 
space dependence of the expectation value (0) and take the wave-function renormalization 
to be constant, = 1. Since the anomalous dimension associated with Z^ is small 
compared to one, see e. g. Ref. [27], our approximation, in which the running of the wave- 
function renormalization is neglected, is well justified for a first study of finite-size scaling. 
The components of the vector (p are labeled according the role the corresponding fields 
are playing in the spontaneously broken regime, (fr^ = (a, tt^, tt^, vr^). We choose the first 
component to be the radial mode in the regime where the ground state of the theory is not 
symmetric under 0(4) transformations: 



The RG flow equation for the effective action according to C. Wetterich [2H] reads: 



where the dimensionless flow variable t is given hj t = ln{k/A) and A denotes a UV cutoff 
at which all couplings are initially specified. The regulator function specifies the details 
of the Wilsonian momentum-shell integrations and has to satisfy certain constraints |2B]- 
Since the choice of the regulator function is at our disposal, we can use it to optimize the 
RG flow [291 ISni [311, |32] . In the following, we employ the optimized regulator function [3T] 



(0) = 0;5^ = (ao, 0,0,0). 



(2) 




(3) 






We then find for the fiow equation for the effective potentials [211 ES] : 
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Here we have replaced the bare masses and couphngs in the inverse two-point functions with 
the scale-dependent quantities. The quantities and are the eigenvalues of the the 
second-derivative matrix of the potential. Note that these quantities still depend on the 
background fields a and vf^. 

In Refs. [Ml EHl [26], the proper-time Renormalizaton Group (PTRG) has been used to 
study finite-volume effects in a quark-meson model. Within the PTRG framework, the RG 
flow equations for the potential can be derived straighfowardly by inserting a cutoff function 
into the Laplace-transform of the one-loop effective action [231 Ell E31 EH ESI ES] • Although 
the PTRG framework does not, in general, yield exact RG flow equations [331 EZ], the PTRG 
cutoff function can be chosen such that the resulting flow equations for the effective potential 
in LPA are identical to Eq. ([s]), which was found in Refs. [311 133] . 

For studying scaling behavior it is convenient to deal with dimensionless quantities rather 
than dimensionful quantities. Therefore we introduce the dimensionless potential m, the 
dimensionless masses m„ and well as the dimensionless field-vector by 



Integrating the flow equation from the UV scale A to A; —> 0, we obtain an effective potential 
in which quantum corrections from all scales have been systematically included. 

Since we are eventually interested in phase transitions in finite volume we need a linear 
term with a source term H in the ansatz for the effective action, which corresponds to an 
external magnetic field in a spin model. In order to solve the RG flow for the effective 
potential U (or m), we expand the potential in local n-point couplings around its minimum 



Uk{a,n^) = ao{k) + a^{k){a^ + vf ^ - (y^{kf) + a2{k){a'' + vf ^ - a^{kff + . . . - Ha , (8) 



where H is the flxed, external symmetry-breaking fleld and all other couplings and the 
minimum are scale-dependent. Since we have absorbed the symmetry-breaking linear term 
into the ansatz for the potential, Uk depends on the flelds cr and vf separately. The condition 



Uk = k ^Uk mlj^ = k "^M^f^ and ipi = k 2 0. . 
Applying these deflnitions to the flow equation (Isl), we obtain 



(6) 





dUk{cr,n^) 
da 







(9) 



(T = (T0(fe),7f2=0 
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ensures that we are expanding around the actual physical minimum. From Eq. (|9]), we find 
that the RG fiow of the coupling ai{k) and the minimum ao{k) are related by the condition 

2ai{k)ao{k) = H . (10) 

This condition keeps the minimum at {a, vr) = {ao{k), 0). The fiow equation of the minimum 
(Jo{k) is thus related to the fiow of the coupling ai{k) in a simple way. 

The RG fiow equations for the couplings aj(fc) can now be obtained straightforwardly 
by expanding the fiow equation ([S]) around the minimum ao{k) and then projecting it onto 
the derivative of the ansatz ([s]) with respect to k. This procedure results in an infinite 
set of fiow equations for the couplings ai{k). In order to solve the set of equations for the 
couplings, we have to truncate our ansatz ([s]) for the potential. In the following, we include 
fiuctuations around the minimum up to eighth order in the fields, i. e. we keep track of 
the running of the couplings ai, 02, as and 04. The resulting finite set of coupled first-order 
differential equations is then solved numerically. From investigations of the convergence 
behavior [271 ES] , we expect that such a truncation is sufficient for our purpose. Below we 
confirm explicitly that the scaling functions for the couplings of interest (ctq and 02) satisfy 
the expected scaling relations, which is strong evidence that the expansion has converged 
sufficiently and the truncation at this order is justified. 

B. Finite volume 

Now we generalize the RG fiow equations in the first part of this section to a finite d = 3 
dimensional Euclidean volume. This is done by replacing the integrals over the momenta in 
the evaluation of the trace in Eq. ([3]) by a sum 

We only consider isotropic volumes, but the approach is not limited to these. Anisotropic 
Euclidean volumes have been used in a study of the quark-meson model in d = 4 dimensions, 
see Ref. [26]. The boundary conditions in the Euclidean time direction are fixed by the 
statistics of the fields, i. e. we must choose periodic boundary conditions in this direction. 
In the present case, there are only three spatial dimensions. In order to be able to compare 
our results with Lattice simulations of 0(A^)-models, we choose periodic boundary conditions 
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in these spatial dimensions. Thus the 3-momenta are discretized as follows: 

4vr2 

P' = ^{nl + nl + nl), (12) 

where rij G Zq for z = 1, 2, 3. The flow equation for the effective potential u can now be de- 
rived straightforwardly. From the flow equation ([s]) together with the regulator function (|4]), 
we obtain 



dtUt 



= -3u + 3Lo{ml,kL) + Lo{ml,kL), (13) 

ERG 



where we introduced the finite-volume threshold functions 

^^(-.-) = i7r^MiEe(i-(¥)V). (14) 



The sum counts the number of lattice nodes located in a three-dimensional ball with radius 
In the limit a; ^ oo (i. e. L ^ oo), the sum can be written as an integral, which yields 
the well-known threshold functions in three dimensions 131): 



' (4vr)ir(| + l) (1 + ^)^+1 ^ ' 

In the preceding part of this section, we have pointed out that a PTRG cutoff function 
can be chosen such that the resulting PTRG flow equations are identical to those obtained 
from the flow equation (|3]) when the optimized regulator Q is employed. Although this 
correspondence is true for inflnite volume, this is not the case for flnite volume. In the 
following, we use the same PTRG cutoff function as in the derivation of the flow equation ([T]) 
in inflnite volume. The flow equation can be derived along the same lines as in Refs. 
and we flnd 



dtUt 



PTRG 



-3m + 3 Qf\ml, kL) + Qf\ml, kL). (16) 



For convenience, we have introduced the (dimensionless) threshold-function 



where dp is the Jacobi-Elliptic-Theta function deflned by 



oo 



dp{x) = ^ e""™' = x-5 +2^x-5e-^. (18) 

71=— oo g=l 
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The first representation in Eq. (18) corresponds to the standard Matsubara summation of 



the momenta. The second representation on the right hand side is obtained by applying 



Poisson's formula to the first representation. Note that the PTRG flow equation (16) is 



identical to the flow equation (13) in the small- volume limit (L — 0), where only the zero 
modes contribute, and by construction also in the infinite- volume limit (L oo). In order to 
obtain the infinite-volume flow equation analytically from the corresponding finite-volume 
flow equation, we use the second representation of the Jacobi-Elliptic-Theta function in 



Eq. (18). This representation actually separates the finite- volume and the infinite- volume 



contributions. When we approximate the Jacobi-Elliptic-Theta function by the first term 



on the right-hand side of Eq. (|18j) and then perform the integration over s in Eq. (17), we 
do indeed recover the flow equation ([T]) for infinite volume. 

On the one hand, physical quantities should not depend on the choice of the regularization 
scheme in the limit k 0. One the other hand, investigating the regulator dependence of 
the results in this limit allows us to check their quality in a particular truncation. In this 



paper, we exploit the difference between the results obtained from Eq. (13) and Eq. (16) in 
order to obtain a theoretical error estimate for our results. In addition, in both schemes the 
restriction to the local potential approximation introduces an additional systematic error. 

III. FINITE-SIZE SCALING 

Critical behavior in the vicinity of a critical point is governed by the presence of long-range 
correlations. A finite volume V = L'^ affects the critical scaling behavior if the volume size 
L becomes comparable to the correlation length ^. According to Fisher's finite-size scaling 
hypothesis [10], observables in the finite- volume system and the infinite- volume system are 
then related by a function that depends only on the ratio of the infinite-volume correlation 
length ^(t, h,L ^ oo) and the linear volume size L. 

We use the customary notation and denote the reduced temperature by t = ^ with 
the critical temperature Tc, and the external symmetry-breaking field hj h = j^. The values 
for Tc and the non-universal normalization constants Tq and Hq which we determined for 
our parameter choices are given in Tab. [1} For completeness, we include the values of the 
critical exponents used in the evaluation in Tab. |TTj For a discussion, we refer the reader to 
Ref. m. 
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A. Finite-size scaling functions 

Applying the scaling hypothesis for example to the order parameter M, the ratio of its 
values in finite and infinite volume is given by 



M{t,h, L oo) V L 

An RG analysis tells us how the couplings t and h have to be changed in order to keep the 
system invariant under a change of the length scale. If the volume is finite in all dimensions, 
the same critical exponents as for infinite volume govern the scaling behavior. Since the 
correlation length behaves as ^ ~ in the absence of the field keeping the ratio of 
correlation length ^ to system size L constant requires to hold tL^^'^ = const. Likewise at 
the critical temperature t = 0, we need to keep the combination hL^^^'^ = const, in order to 
preserve the ratio of correlation length and system size. 

This completely specifies the behavior of thermodynamic observables in the vicinity of the 
critical point in a finite volume system, apart from possible scaling corrections. For example, 
the order parameter M{t, h, L) as a function of temperature t, external symmetry-breaking 
field h and volume size is expected to behave as [IQ] 

1 



M{t,h,L) = L-^'" 



(20) 



where u is the critical exponent associated with the first irrelevant operator in the vicinity 
of the critical point. The additional terms spoil the finite-size scaling behavior for small 
volume size. They need to be accounted for explicitly and removed to isolate the universal 
finite-size scaling function. 

In general, the finite-size scaled order parameter L^^^M{t^ h, L) is a function of the finite- 
size scaled variables tL^^'^ and hL^^/'^ . In order to analyze the results efficiently in terms of 
only a single variable, we parametrize the scaling functions in term of the scaling variable 
z = This accounts for the usual critical scaling behavior in infinite volume. We 

keep the value of z fixed and use the finite-size scaling variable 

h = /iL^'^'^)/^ (21) 



to parametrize the finite-size scaling behavior. The scaling variable h is not dimensionless 
but retains an explicit scale dependence through the length scale L. While we have not done 
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so in the present analysis, it is possible to remove this explicit non-universal scale as well: 
Using the correlation length ^ as the natural length scale, we can form the dimensionless, 
universal combination 

^ = M,(,/..L)L, (22) 

which still depends explicitly on t, h (and L). This relation can be used to identify a unique 
length scale by considering the limit L — > oo and /i — > 0. In this limit, the correlation length 
behaves according to 

M,{t) = Cotr (23) 

The non-universal normalization constant Co has the dimension of an inverse length scale 
and plays the same role for the " coupling" L as the constants Tq and Hq do for the couplings 
T and H. Using this constant, the universal dimensionless combination M^t^L — > oo)L 

becomes 

where £ = CqL is the dimensionless length required for a truly universal description of the 
finite-size scahng behavior. With our choice of parameters for a UV cutoff scale A = 1 GeV, 
we find Co = 2.393(2) MeV, corresponding to a length scale 

= 82.510(65) fm. (25) 

Co 

We stress that all explicitely given length scales must be normalized with this factor for a 
completely scale-independent comparison. 

In the following part of this section we review some properties of the finite-size scaling 
functions and make contact with the scaling results in infinite volume 

Provided the finite-volume corrections can be neglected or removed, the finite-size scaled 
order parameter is given by the universal finite-size scaling function 

L^/"M(t, h, L) = Qm{z, /iL(^^)/"), (26) 

where the scaling variable z is used to parameterize the dependence on the temperature t. 
For large volumes, the order parameter must asymptotically approach the infinite- volume 
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scaling limit 

Urn M{t,h,L) = h^/^f{z). (27) 

This dictates the behavior of the finite-size scaling function Qm{z, /iL^^^^/^) for large values 
of its second argument, hL^l^^^^'^. In order to recover the infinite-volume scaling law, the 
magnetization must behave as 



lim M{t,h,L) = lim L-^/"Qm(^, /ii^^^^^^") = lim L"^'"' ihL^I^^^/") ' f{z) 

'y—>oo L—>oo L^oo 

= h'/'f{z). (28) 



Therefore we expect to find 

QM{z,h) - f{z)¥/^ for h^oo (29) 

for the scaling function. While this is a very useful check for the consistency of the finite- 
size scaling results with the infinite-volume scaling behavior, it is only of limited practical 
interest for a finite-size scaling analysis: In the region where the finite-size scaling function 
displays this behavior, the actual results for the order parameter have already converged to 
the infinite-volume result, and there are no large observable finite-size effects. In this region 
the results for the order parameter from different volume sizes coincide, apart from possible 
exponentially small corrections, but the finite-size scaled results L^/'^M(t, /i, L) no longer 
do, see e.g. Fig. [Tjfor 2; = and Fig. [2]for z = Zp. Recall that z = means that we are 
sitting at the critical temperature whereas z = Zp means that we are sitting at the peak 
value of the susceptibility. 

Since the longitudinal susceptibility x is given by the derivative of the order parameter 
with respect to the external symmetry-breaking field if, the infinite- volume scaling functions 
f{z) for the order parameter and /^(-z) for the susceptibility are related. In infinite volume, 
one finds the relationship 



I3h 

As in infinite volume, the finite-size scaling functions for the order parameter and for the 
susceptibility are related to each other: Taking the derivative of the finite-size scaled order 
parameter leads directly to the corresponding expression for the susceptibility. The finite- 
size scaling function Qm{z, hL^^^") depends on two variables z and h. Taking the derivative 
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with respect to /i, derivatives with respect to both variables appear, which we denote by 

QM°\z,h) ^ -^QM{z,h) and QM^\z,h) ^ -^QM{z,h). (31) 
It is further useful to remember that 

dh ^ dh (55 h 

We obtain the finite-size scaling function for the susceptibility by taking the derivative 

of the finite-size scaling function for the order parameter: 



1 jpS/u 



(32) 



(35 hL^^I^ 

Using the scaling relation 7 = (5{5 — 1), one finds that the finite-size scaled susceptibility is 
given by the scaling function 

HoL-'^/^x{t,h,L) = Q'-^'\z,hL^'/n - j^j-^Q^i^'\z,hL^'/n ^ Q^{z,hL^'/'^). {33) 

Since we determine the scaling functions for fixed values of z as a function of the finite-size 
scaling variable hL^^^^'^, we can test this relationship between the finite-size scaling functions 
most easily for the case z — 0, where the second term vanishes. Here we assume that the 
derivative of Qm with regard to z is bounded for 2; —> 0, or that HuIz-^qz Q^^^\z,h) — 0, 
which is indeed the case in our results. 

For large volumes, once again the finite-size scaling function must coincide with the 
scaling functions in infinite volume, 

hm xit, h, L) = hm L^/-' [hL^"''f"~' ^f^{z) = h'/'-'^f^iz), (34) 

L— >0O i— >0O IlQ riQ 

where we used the scaling relation (3{5 — 1) = ^ again. Consequently the scaling function 
must behave asymptotically for large values of h as 

Q^iz^h) = h'/'-'^f^{z), (35) 

where the scaling function fx{z) and the normalization constant Hq are known from the 
infinite- volume result. 
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B. Scaling behavior as a function of the correlation length 



The principal problem in comparing universal scaling functions to actual data is the 
determination of the non-universal normalization constants. The determination of these 
scales is essential for a comparison of the scale-free universal results, e. g. the constants in 
a Lattice QCD simulations are in general different from those in a lattice simulation of an 
0(4) spin model. But the determination from a limited data set is difficult, since such a set 
is often obtained in the presence of strong symmetry-breaking fields and in small volumes. 
These may lead to non-universal scaling corrections where details of the short-range physics 
enter. 

The most useful results are thus those that rely only on the long-range properties of the 
system in question. Dimensionless variables into which only the actual observables enter and 
for which any normalization with regard to short-range physics is unnecessary are actually 
the best candidates for a practical evaluation. 

Considering the discussion of finite-size scaling above, the most natural variable for plot- 
ting finite- volume results is the ratio of the infinte- volume correlation length to the system 
size 



But the infinite-volume correlation length is itself not directly measurable either, so it is 
necessary to use the finite-volume correlation length, which can be measured in the finite- 
volume system. 

Since the correlation length is an observable, we can apply the finite-size scaling analysis 
to the correlation length itself. According to the general scaling analysis, for any observable 
O, we have up to scaling corrections 



where k is the critical exponent associated with the operator O. In particular, we have 
K = V for the correlation length itself and the relation then becomes 



Assuming the scaling function in this relationship can be inverted, we can express the ratio of 
the infinite-volume correlation length and the volume size as a function of the finite-volume 



^{t,h, oo) 



L 




(36) 




(37) 
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correlation length and the volume size: 



^{t, h, oo) _ / ^{t,h,Ly j ^gg^ 



L \ L 



Substituting this relation into the scaling hypothesis Eq. (19), the scaling functions can be 
expressed as a function of the argument ^^^''^'^^ involving the finite- volume correlation length. 
This form is very useful for actual comparisons of different physical systems. 



In Sect. VA, we will show results for the finite-size scaling function Qm{^/L) for the 
order parameter at selected values of the scaling variable z {z = and z = Zp = 1.3155) 
in this form. Since the correlation length and the susceptibility are trivially related in the 
present approximation, x = ^/M^ = no additional information is gained by determining 
the finite-size scaling function of the susceptibility as a function of the ratio C,/L. 

A different comparison scheme which also eliminates the need for additional normalization 
constants is based on the comparison of the system for two different volume sizes, L and sL 
with s G M. This approach is widely used in condensed matter physics, see e.g. Refs. [H], 
2l 03]. The most common choice for the volume ratio is s = 2. 
Starting from the inverted scaling relation for the finite-volume correlation length 



Eq. (38), one observes that for a volume of size sL 



^{t, h,oo) ^ h,oo) ^ l^_i h,L) \ ^gg^ 

sL s L s ^ \ L 



Inserting this expression into the scaling relation Eq. (36) and forming the ratio for system 
sizes L and sL, one obtains 

r (1 f-l ( (it,h,L)\\ 

oit, K sL) _ ^^/ Joy-sh y—)) ^_ i{t,KL) \ _ ^^^^ 



Scaling corrections to this relation are of the order and L~'^, where uj is the exponent 
associated with the first irrelevant operator. The ratio of the values of an observable for 
two systems sizes L and sL therefore defines a new scaling function Fo which depends 
on the ratio s and the dimensionless correlation length ^{t, h, L)/ L. This scaling function 
can be used for direct comparisons without any need for an additional determination of 
normalization constants. 



In Sect. |VC[ we discuss the results for this finite-size scaling function for the order 
parameter and for the correlation length, which corresponds to the susceptibility, for selected 
values of z {z = and z = Zp = 1.3155). 
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IV. DETERMINATION OF THE FINITE-SIZE SCALING FUNCTIONS 

Ideally, the results for different values of the volume size in the finite-size scaling region 
would coincide perfectly after rescaling. But as we have already observed for infinite- volume 
scaling [23], the external symmetry-breaking field h has to remain small in order to keep 
corrections to scaling small. This means in particular that the masses of the critical fluctu- 
ations must remain very small compared to the UV cutoff scale A. On the other hand, in 
order to fully explore the finite-size scaling region, we must decrease the correlation length 
sufficiently to restore the infinite-volume results, which implies an increase in the masses of 
the fluctuations. 

Therefore corrections to the finite-size scaling behavior are already large for results from 
relatively large volumes, compared to the scales set by our choice of parameters. In order 
to have a sufficiently large set of results to extract the scaling behavior, we also take results 
with non-negligible scaling corrections into account. 

As outlined above, including the first explicitly volume- dependent correction term, the 
order parameter can be expanded as 

1 



M{z,h,L) = L-l^'' 



(41) 



In order to extract the scaling functions for the order parameter, we start from the assump- 
tion that the deviation from the perfect scaling behavior contained in the scaling function 
QM{z,h) is completely determined by the function Q^^^ in Eq. (41). We will confirm that 
this assumption is sufficiently well satisfied by our results. 

We can remove the leading-order scaling function by comparing the rescaled order param- 
eter for different volume size, but at the same value of the scaling variables z and h = hL^^^'^ . 
This determines the coefficient function Q^^^\z, h) of the non-universal correction term. For 
two different volume sizes Li and L2, we choose values for the external fields hi and /12 such 
that 

h* = hi = hiLf^^" = h2 = /isi^f (42) 
Since the finite-size scaling function depends only on z and h, the universal term Qm{z, h) 
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drops out in the difference of the two values and we are left with 

= Q^^{zrh*)^{s--l) + ... (43) 

with s = L2/L1. By fitting this expression as a function of the size ratio s to results for 
different volume sizes, the coefficient function as well as the exponent u can be determined. 

Scaling corrections due to large values of the external symmetry-breaking field and the 
small volume are both present in our results. Since these scaling corrections become pro- 
gressively larger the further a way one moves from the critical temperature and the smaller 
the volume, we have restricted the results for large values of the variable z to larger volume 
sizes. For \z\ < 1.0, we take volume sizes from L = 100 fm down to L = 1 fm into account, 
while for \z\ > 1.0 we only use the volumes from L = 100 fm down to L = 10 fm to extract 
the finite-size scaling functions. 

V. RESULTS AND DISCUSSION 

In this section, we present the results from our investigation of the finite-size scaling 
behavior of the 0(4) model in c? = 3 dimensions and from our determination of the finite-size 
scaling functions. We start with a discussion of the finite-size scaling behavior of the order 
parameter, illustrate the extraction of the scaling functions, and propose a parameterization 
for these functions for the limits of small and large values of the scaling variable h. We 
continue with an analogous discussion of the finite-size scaling behavior of the susceptibility. 
We finally show results for the finite-size scaling functions QM{z,h) and Q^{z,h) for a 
selected set of values for the scaling variable z as a function of h and discuss effects of 
our approximation scheme. In the last part of this section, we briefly discuss the scaling 
functions as functions of the dimensionless correlation length ^/L. 

A. Finite-size scaling function for the order parameter 

The behavior of the order parameter as a function of the external symmetry-breaking 
field h for different volume sizes is shown in the first panels of Fig. [T]for 2; = 0, i. e. at the 
critical temperature, and Fig. [2]for z = Zp, i. e. at the peak value of the susceptibility. 
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Figure 1: Finite-size scaling behavior of the order parameter M at the critical temperature, z = 0. 
Shown is the unsealed result for M as a function of h for different volume sizes (first panel) , and the 
finite-size scaled order parameter L^^'^M as a function of hL^^^'^ for the same values of the volume 
size (second panel). For large values of h, where the correlation length is small, the unsealed results 
for different volume sizes all converge towards the same infinite- volume limit (first panel) . For small 
volume sizes, the corrections to the ideal scaling behavior become considerable. The black dashed 
line (second panel) is the result for the scaling function obtained from all results by removing the 
scaling corrections. 
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Figure 2: Finite-size scaling behavior of the order parameter M at the peak of the longitudinal 
susceptibility, z = Zp = 1.3155. Shown is the unsealed result for M as a function of h for different 
volume sizes (first panel), and the finite-size scaled order parameter L^^'^M as a function of hL^^/'^ 
for the same values of the volume size (second panel) . The black dashed line (second panel) is the 
result for the scaling function obtained from all results by removing the scaling corrections. 
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Two separate regimes can be distinguished in the double-logarithmic plots: For large 
values of h, the masses of the fluctuations are large and hence the correlation length is much 
smaller than the volume size. Even the results for small volume size approach asymptotically 
the infinite-volume limit. As expected from the infinite-volume scaling behavior, the order 
parameter behaves as ~ h^^^ and the slope of the curves for large fields is determined by the 
value of the critical exponent 6. For small values of h, where the correlation length becomes 
large, the deviations from the infinite-volume behavior are equally clear: For small volumes, 
a small correlation length can be of the order of the volume size and the deviations occur 
already at large values of h, whereas for large volumes the correlation length needs to be 
large for significant finite-size effects which appear only for small values of h. 

In the second panels in both Figs, [l] and [2| the rescaled order parameter ML^^'^ is 
plotted as a function of the scaling variable hL^^^^. In the finite-size scaling region, the 
curves for different values of the volume size collapse almost perfectly onto a single curve. 
The agreement becomes worse with decreasing volume size, which can be explained by 
the presence of scaling corrections. When we assume that these scaling corrections can 
be described by the corrections due to only the first irrelevant operator, we can use the 



expression from Eq. (43) to fit the results. After subtracting the corrections, we obtain a 
result for the universal leading-order scaling function QM{z,h) (shown as a black dashed 
line in the figures). 

From an RG determination at the fixed point [321 SH SS] we expect for the exponent u 
associated with the corrections the value 

CO = 0.7338. (44) 

For z = 0, we find from a fit to our rescaled results for L = 10 fm to L = 100 fm for the 
order parameter in the finite-size scaling region with h < 1.0 the value 

CO = 0.7443(300), (45) 

which is compatible with the fixed-point determination. Since we fit to a numerical result 
that also includes contributions from additional irrelevant operators, our result for co should 
be regarded as an effective value. Having established the consistency of the scaling correc- 
tions in our results with the RG prediction, we use in the following the more accurate value 



(44) for the evaluation. 
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We have determined the scahng functions Qm{z, hL^^^'^) for a range of z-values. Plots 
for selected values are shown below in Fig. [3]for z > 0, and in Fig. |4]for z < 0. For values 
\z\ < 1.0 we have used results from volume sizes L = 1 fm to L = 10 fm in steps of 1 fm and 
from L = 10 to L = 100 fm in steps of 10 fm. For values \z\ > 1.0, we have only included 
results from volume sizes from L = 10 to L = 100 fm in steps of 10 fm. For small volumes, 
the value of h has to be increased too much in order to access the finite-size scaling region 
and scaling corrections become too large. 

We have not found a global parameterization for the scaling function, but we find that 
for small values of hL^^^^ our results are fitted well by the parameterization 

q£)(z, hL"'/") = c{z){hL^'/''y'^'\ (46) 

Both the coefficient c{z) and the exponent t{z) vary with z. The values we determined for 



the fitting parameters can be found in Tab. |III| and Tab. IV for \z\ < 1.0 and \z\ > 1.0, 
respectively. We include in the tables the approximate values for hL^^^'^ below which the fit 
is applicable. For comparison to the actual scaling functions, we also show the fits in the 
first panels of the Figs. [3] and |4} 

For large values of /iL^*^/^, the parameterization is constrained by the asymptotic behavior 
of the scaling functions for infinite volume: 

qSH^, hL^'/^ = c^{z){hL^"y-. (47) 

We expect the exponent Too to be independent of z and to have the value Too = l/<5 = 
0.2011(1) for our values for the critical exponents. The coefficient Coo{z) retains its z- 
dependence and is expected to coincide with the value of the infinite-volume scaling function, 
Coo{z) = f{z). The results from the fits are summarized in Tab. |V]for \z\ < 1.0 and Tab. VI 
for \z\ > 1.0. For comparison with c{z) and Too, the corresponding values for f{z) are 
included in the tables. The values for the coefficient Coo{z) and f{z) agree within 1%, 
whereas the agreement of the exponent Tqo with the expected value is only fair. 
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Figure 3: Finite-size scaling functions for the order parameter (left column) and the susceptibility 
(right column) as a function of hL^^/'^ for fixed z > from PTRG (red solid line) and ERG with 
optimized cutoff (blue dashed line). For z = 0.1, the results of the fits to the asymptotic behavior 
for small and large values of the scaling variable are included (black dashed lines). The asymptotic 
behavior is described well, but the parameterization is clearly not applicable in the intermediate 
region. 
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Figure 4: Finite-size scaling functions for the order parameter (left column) and the susceptibility 
(right column) as a function of the finite- size scaling variable hL^^/'^ for fixed z < 0. Shown are 
results obtained from both the PTRG (red solid lines) and the ERG with optimized cutoff (blue 
dashed lines). Both cutoff schemes agree well in the finite-size scaling region {hLl^^l'' < 1), and 
deviations appear only for larger values, where the functions are already determined by the infinite- 
volume behavior. The black dashed lines included for z = represent the fits to the asymptotic 



23 



B. Finite-size scaling function for the susceptibility 

The scaling behavior of the susceptibihty x as a function of the field h is shown in Figs. [5] 
and |6] for z = and z = Zp = 1.3155, respectively. We can once again distinguish an 
asymptotic region for large h, where the correlation length is much smaller than the volume 
size and the results converge towards the infinite- volume limit, behaving as ~ 

For 

small values of h, the finite-size effects are strong and the susceptibility tends towards a 
constant value which depends on the volume and varies by about two orders of magnitude 
when the volume is changed by a factor of 10. 

In both Figs, [s] and [g] the results for the rescaled susceptibility L~'^'''^x a^re shown in the 
second panels as a function of the scaling variable hL^^^'^. As for the order parameter, the 
collapse of the rescaled results onto the scaling function is quite good, although not perfect 
due to scaling corrections for smaller volume size. We use once again the RG description of 
these scaling corrections to extract the leading-order scaling function, which is depicted by 
a dashed black line in the figures. 

Scaling corrections from the results for the rescaled susceptibility for z = at h < 1.0 
lead to a value of the exponent 

uj = 0.7384(400) (48) 



which is in agreement with the result Eq. (45) from the order parameter, and with the result 



Eq. (44) from the fixed point determination [33l HH HS]. As for the order parameter, we use 



the more accurate result (44) for all further evaluations. 

The determination of the finite-size scaling function and the values included in the eval- 
uation are the same as outlined above for the order parameter. 

As for the order parameter, we have not found a satisfactory global parameterization for 
the scaling function Q^i^, h). However, we are able to parameterize the functions for small 
and large values of h separately. 

For small values of /i, we find that the parameterization d{z) + c{z){hL^^^'^y works 
sufficiently well in the range of 2;- values that we consider. The parameterization of the scaling 
function for the susceptibility is in principle constrained by that of the order parameter. 
Because the exponent r of the parameterization for the order parameter is close to one, 
r— 1 ^ 1, the leading term of the derivative of the scaling function Qm{z, h) with respect to h 
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is approximately constant, which is consistent with the ansatz for the parameterization that 
we have chosen here for the susceptibihty. Beyond that, the compatibihty of the results for 
susceptibility and order parameter cannot be checked very well using this parameterization. 
We will discuss this in more detail below and demonstrate for the special case z = that 
our results for the order parameter and the susceptibility are indeed consistent. The values 



for the fit parameters determined for small values of h are given in Tab. VII for small values 



of 1^1 < 1.0 and in Tab. VIII for values of \z\ > 1.0. The last column in the table lists 
the approximate values below which the parameterization is valid. For comparison to the 
scaling functions, the fits are shown in the first panels of Figs. |3]and|4j 

For large values of the scaling variable, the parameterization of the scaling function is 
once again constrained by the known behavior in the infinite- volume limit. We use the 
parameterization Coo{z){hL^^^^Y°° for large values of h. We expect that the exponent is 
independent of z and takes the value Too = 1/5 — 1 = —0.7989(1) with our value for S. 
For the coefficient, we expect that it coincides with the infinite- volume scaling function for 
the susceptibility, Coo{z) = -^^f^{z). The results of the fits are summarized in Tabs. 
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and X for \z\ < 1.0 and \z\ > 1.0, respectively. For comparison, values for HqCoo{z) and 
fx{z) are included in the same tables. The agreement for the coefficient Coo{z) with the 
scaling function /^(-z) is not as good as the agreement in the corresponding relation for the 
order parameter, although still quite good close to the critical temperature {z = 0) and 
the susceptibility peak {z = Zp = 1.3155). The agreement of the exponent Tqc with the 
expectation is good for small values of z, but for both parameters the agreement becomes 
progressively worse when one moves to large absolute values of z. 



As discussed in Sec. Ill, the finite-size scaling functions Quiz^h) and Qx{z,h) of the 
order parameter and the susceptibility are connected, and in fact the scaling function for 
the susceptibility should be completely specified by the one for the order parameter. Due 
to the dependence of the fit parameters on the scaling variable z and the limited number of 
2;-values for which we determined these parameters, we cannot test this relation in terms of 
the parameterization. In addition, because of the additional 2;-dependence we also cannot 
test the relation numerically for arbitrary value for z. However, as we have seen from 



Eq. (33), for the special case z = additional 2;-dependent terms drop out and we can check 
the validity of this relation. The comparison between the actual scaling function Q^{z,h) 
obtained from the susceptibility and the derivative of the scahng function Qm{z, h) is shown 
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in Fig. [7] as obtained from our PTRG approach. We find that the scahng function obtained 
from the susceptibihty (sohd red hne) and the derivative of the scahng function for the 
order parameter (dashed black hne) agree very weh, and are indistinguishable on the scale 
of the plot. This demonstrates that the RG flow equations capture the scaling behavior 
consistently and that the results do indeed satisfy the expected scaling relations, as we have 
already seen in the application to infinite volume [23] . With regard to the truncation of the 
RG flow, the agreement demonstrates that the truncation of the RG flow is sufficiently large 
for our calculation. 

The results for the scaling functions for both the order parameter and the susceptibility 
shown in Figs. [3] and |4] have been obtained from an RG calculation using a proper-times 
cutoff (PTRG) (red lines) and from an ERG calculation using an optimized cutoff (blue 
lines). For infinite volume, both calculations coincide and the results are exactly the same. 
In addition to the agreement in infinite volume, the result again coincide for very small 
values, where all fluctuations are cut off by the finite volume size. In finite volume, the 
threshold functions that enter into both calculations differ, see the discussion in Sec. |Tl} 
Due to truncation effects, which affect the finite-volume RG fiows in both calculations in 
slightly different ways, the results do not coincide in the intermediate-volume region (where 
2.5 < L/^ < 5). Note that the double-logarithmic scale in Figs. [3] and |4] magnifies the 
apparent difference between the results for large volume. For asymptotically large external 
symmetry breaking fields (mean-field limit), we find agreement of the results from both 
approaches again. Overall, we observe that the agreement is better for the order parameter 
than for the susceptibility, and we conclude that the truncation errors affect the higher-order 
couplings that enter into e.g. the susceptibility more strongly than the order parameter. The 
difference between the methods can be used to assess the truncation errors and provides a 
theoretical error estimate. 
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Figure 5: Finite-size scaling behavior of the longitudinal susceptibility at the critical temperature, 
z = 0. The susceptibility for different volume sizes is shown as a function of h (first panel). For 
large values of h, the results for different volume sizes all converge towards the infinite-volume 
limit. The finite-size scaled susceptibility L^'^^'^x for the same volume sizes is shown as a function 
of hL^^/^ (second panel). For small volume sizes, the corrections to the ideal scaling behavior are 
of considerable size. The result for the scaling function obtained from all results by removing the 
scaling corrections is shown as a black dashed line. 
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Figure 6: Finite-size scaled longitudinal susceptibility at the peak value, z = Zp = 1.3155. The 
susceptibility x for different volume sizes is shown as a function of h (first panel) . For large values 
of h, where the correlation length is small, the results for different volume sizes all converge towards 
the same infinite- volume limit. The finite-size scaled susceptibility L~^l^x is shown as a function 
of hL^^I^ (second panel). Corrections to the scaling behavior can be seen for the small volume 
sizes. The scaling function obtained by removing the scaling corrections is represented by the black 
dashed line. 
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Figure 7: Comparison of the finite-size scaling function Q^{z, h) determined from the susceptibihty 
(red sohd hne) and Q^^^^\z, Ji), the derivative of the scahng function for the order parameter with 
respect to h (black dashed line), for z = 0, as obtained from the PTRG. As expected, they agree 
well and are almost indistinguishable on the scale of the figure. This provides strong evidence that 
the truncation of the RG flow is sufficient to account for the relevant scaling behavior. 
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C. Scaling results as a function of the correlation length 



As discussed in Sec. Ill , for comparisons of different systems in the same universality class 
it is advantageous to use dimensionless ratios since it eliminates the need for a determination 
of non-universal normalization constants. For this reason, we briefly present results for the 
scaling behavior as a function of the dimensionless finite- volume correlation length ^(L)/L 
in this section. 

Using the relation x{^jh.,L) = C,'^{z,h,L) for = 0, we translate the scaling function 
for the order parameter Qm{z, h) into a function Qm{z,^{L) / L) of ^{L)/L. The results for 
z = Q and shown in Fig. pj Since the susceptibility is bounded for a given volume, 

the correlation length is bounded as well, and we observe that this bound is approximately 
at ^{L)/L ^ 0.5. 

Alternatively, one can also calculate the ratio of an observable for two different volume 
sizes L and sL, which differ by a fixed factor s, and plot the results as a function of the 
dimensionless correlation length ^{L)/L for the smaller volume. The results of this procedure 
for the order parameter M{2L)/M{L) and for the correlation length itself ^{2L)/^{L) are 



shown in Figs. |9] and 10 We have again chosen the values 2 = and z = Zp. The plots 
are generated by keeping the scaling variable z fixed while varying h, which changes the 
correlation length. As expected, we find that for both observables the ratio 0{2L)/0{L) 
tends to one in the infinite- volume limit ^{L)/L 0, where the correlation length for a 
given external symmetry breaking field becomes much smaller than the volume size and 
the finite-size effects disappear. In small volumes, corrections to the scaling behavior are 
unfortunately large, and a systematic extrapolation is more difficult here, since the scaling 
functions are given in parametric form. But it is clear that the results converge to a universal 
result for these scaling functions for large volume size, as well. 

In particular for the correlation length, there is an additional check that helps us to 
evaluate the convergence behavior: As we have already noticed in the discussion of Fig. |8| 
the correlation length is bounded by the volume size (at apparently ~ 0.5L). For growing 
correlation length, the ratio ^(2L)/^(L) in volume sizes 2L and L must eventually converge 
to ^(2L)/^(L) = 2, and with increasing volume size, we observe exactly this convergence to 
the expected value. 

We also observe that the shapes of the scaling functions for the ratios vary with the value 
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Figure 8: Finite-size scaling function QmH/L) for the order parameter as a function of the dimen- 
sionless finite-volume correlation length, ^/L. The results shown are for the critical temperature 
{z = 0, red solid line) and for the peak of the susceptibility {z = Zp = 1.3155, blue dashed line). 
Plotting the scaling function as a function of the dimensionless quantity i/L, we eliminate the need 
for the determination of an additional non-universal length normalization factor. We find that the 
finite-volume correlation length is bounded at approximately ^/L ~ 0.5. 

of the scaling variable which could also conceivably help to distinguish between different 
universality classes. 



31 




Figure 9: Ratio of the order parameter M{2L) /M{L) for systems with the volume sizes 2L and L, 
shown as a function of the ratio of the finite- volume correlation length and the system size, ^{L)/L. 
Here the results are given for z = (at the infinite-volume value for the critical temperature) (first 
panel) and z = Zp = 1.3155 (at the peak value of the susceptibility) (second panel). Corrections 
to scaling are large for small volume sizes (1 — 2 fm for our choice of parameters), but the results 
converge for large volume size (~ 100 fm). 
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Figure 10: Ratio of the finite-volume correlation lengths ^{2L)/^{L) for systems of the volume 
sizes 2L and L als a function of the ratio ^{L) / L of the finite- volume correlation length and the 
system size for the smaller system. Since the finite-volume correlation length cannot exceed the 
system size, for the given choice of ^{2L)/^{L) the ratio is bounded by 2 from above. For large 
volume size {S,{L)/L 0), the ratio must tend to unity. 
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VI. CONCLUSIONS 

Using functional RG methods, we have calculated the finite-size scaling functions for the 
0(4) model in three dimensions for a wide range of temperature, symmetry breaking and 
volume size. The behavior of the scaling functions allows us to identify two distinct regions: 
for large values of the correlation length and small values of the finite-size scaling variable, 
finite-size effects are large, and the behavior of the observables is described by finite-size 
scaling. Above a certain value for the scaling variable, where the correlation length is 
sufficiently small, the results converge rapidly towards the infinite-volume behavior, and 
the scaling functions exhibit asymptotic behavior consistent with predictions from infinite- 
volume scaling. 

We have expanded the effective action in terms of local interactions. As discussed in 
|23j . this introduces a systematic error into the values for all critical exponents used in the 
evaluation of the results. We have used two different RG schemes in our calculation in order 
to assess truncation effects on our results, which allows us to estimate a theoretical error. 
In the infinite- volume limit, the results for both calculations coincide, which is also the case 
for very small volumes, where all fiuctuations are cut off. 

Deviations from the scaling behavior in our results can be understood by RG arguments as 
scaling corrections, and we have checked that the behavior is consistent with that predicted 
by the RG. We have removed these corrections, which appear primarily for very small 
volumes and very large values of the external symmetry-breaking fields, and extrapolated 
from our results to the universal finite-size scaling functions. We provide parameterizations 
for our scaling functions separately for the finite-size scaling region at small values of the 
scaling variables and for the asymptotic region at large values of the scaling variables. 

For selected values of the infinite-volume scaling variable 2;, we have confirmed that the 
expected scaling relations between the scaling functions for the order parameter and the 
susceptibility hold for our results. We have also briefiy presented some results for scaling 
functions in terms of dimensionless ratios of observables, as a function of the dimensionless 
correlation length which allow direct comparisons of different systems in the same 

universality class without any need for additional normalization. 

Our results are useful for a comparison to Lattice QCD simulations. In order to perform 
such a comparison for the scaling of the order parameter (or the susceptibility) to results 
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for the chiral condensate (or the chiral susceptibihty) from Lattice QCD as a function of 
the volume size, a determination of the proper normahzation constants for the reduced 
temperature, the external symmetry breaking field and the correlation length for the QCD 
results is necessary. Once these constants are known, a direct comparison is possible. This 
will help to shed more light on the question whether the QCD phase transition falls into the 
0(4) universality class or not. Such a detailed comparison is work in progress \19\. In this 
respect, we would like to remark that our results at the present stage appear to suggest that 
- at least with regard to a finite-size scaling analysis - current Lattice QCD simulations 
might be in the asymptotic region, where the correlation length is much smaller than the 
lattice size and where large scaling effects cannot be observed. This is consistent with the 
observations previously made in [15]. For a meaningful scaling analysis of Lattice QCD 
results in terms of c? = 3 universality classes, the correlation length ^ on the lattice has to be 
much larger than the length scale set by the inverse temperature. At the same time, it must 
be of the order of the spatial extension L for a finite-size scaling analysis. If these conditions 
cannot be satisfied simultaneously, then the universality class of the QCD phase transition 
cannot be determined from the comparison of Lattice QCD results to scaling behavior as it 
is observed in finite-volume d = 3 models. 
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Appendix A: TABLES 

Table I: Values of the non-universal normalization constants for our parameter set used in the 
evaluation. Following the usual convention, they are defined from the normalization conditions 
M{t,h = 0) = {-tf = (j^y for /i = and t < 0, and M{t = 0,h) = h^/^ = [jr^^'^ for t = 0. 
For details, see [23]. Additionally, we use here £,{t,h = 0) = ^t~'^. the overall scale is set by our 
choice A = 10^ MeV. 

Tc To Hq Cq 

13.682 368165 072 75 MeV 0.014916(5) MeV 6.032(10) MeV^/^ 2.393(2) MeV 



Table II: Values of the critical exponents used in the evaluation. We use the values determined 
in [25j in local potential approximation. See Ref. 123^ for a discussion of the accuracy and effects 
of the approximation. 

P V 6 If] 



0.4030(30) 0.8053(50) 4.9727(5) 1.606(10) 
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Table III: Parameterization of the order parameter finite-size scaling function QM{z,hLf^^^'^) for 
small values of the scaling variable hL^^/^ , determined from the results from L = 1 — 10 fm and 
from L = 10 — 100 fm. The last column indicates the value up to which this parameterization is 
approximately valid. 

z c{z) t{z) {hLf^^/") 

-0.5 1.708(39) 0.9578(144) 0.58 
-0.2 1.585(32) 0.9536(115) 0.57 
-0.1 1.524(25) 0.9455(116) 0.57 

0.0 1.464(31) 0.9371(118) 0.58 

0.1 1.405(30) 0.9284(121) 0.58 

0.2 1.347(30) 0.9194(125) 0.58 

0.5 1.140(28) 0.8734(152) 0.61 



Table IV: Parameterization of the order parameter finite-size scaling function QM{z,hL^^/^) for 
small values of the scaling variable hL^^^'^, determined from the results from L = 10 — 100 fm. 
The last column indicates the approximate value of the scaling variable hL^^^'^ below which the 
parameterization is valid. 

z c{z) t{z) {hL^^I") 

-5.0 4.472(187) 1.1067(227) 0.42 
-4.0 3.953(140) 1.0991(191) 0.43 
-3.0 3.3848(983) 1.0848(158) 0.45 
-2.0 2.6825(726) 1.0444(156) 0.50 
-1.0 2.0647(440) 1.0043(122) 0.54 

0.0 1.4396(308) 0.9340(121) 0.59 

1.0 0.9012(289) 0.8302(179) 0.56 

1.3155 0.7890(315) 0.8069(207) 0.50 

2.0 0.6453(499) 0.7887(314) 0.31 
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Table V: Parameterization QM{z,hLf^^^^) = Coo{z){hL^^/'^Y'^ of the order parameter finite-size 
scaling function, for large values of the scaling variable where the behavior becomes asymptotic, 
determined from the results from L = 1 — 10 fm and from L = 10 — 100 fm. Prom the asymptotic 
properties of the scaling function, one expects that Coo(^) = f{z) and that t^o = ^/5\s independent 
of z. For comparison, the value f{z) of the infinite-volume scaling function and 1/5 are also listed. 

Z Coo(z) f{z) Too 1/5 

-0.5 1.14303(39) 1.15371 0.22611(15) 
-0.2 1.05861(17) 1.06353 0.22807(7) 
-0.1 1.02603(26) 1.03179 0.23014(11) 

0.0 0.99587(25) 0.99916 0.23113(10) 0.2011(1) 

0.1 0.96286(33) 0.96564 0.23343(14) 

0.2 0.93380(28) 0.93122 0.23337(13) 

0.5 0.83478(19) 0.82288 0.23907(10) 
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Table VI: Parameterization QM{z,hLf^^^'^) = Coo{z){hLf^^^'^Y°° of the order parameter finite-size 
scaling function for large values of the scaling variable hL^^^/'^ , determined from the results from 
L = 10 — 100 fm. Prom the asymptotic properties of the scaling function, one expects that 
Coo (-2) = f{z) and that Too = 1/5 is independent from z. Por comparison, f{z) and 1/5 are also 
listed. 



z 


^\jO \ ^ J 


fiz) Tr^ 


1/5 


-5.0 


2.0303(8) 


2.0201 0.2084(2) 




-4.0 


1.8745(7) 


1.8734 0.2101(2) 




-3.0 


1.7013(6) 


1.70843 0.2116(2) 




-2.0 


1.5028(5) 


1.51793 0.2131(1) 




-1.0 


1.2705(3) 


1.28915 0.2145(1) 


0.2011(1) 


0.0 


0.9828(1) 


0.99916 0.2205(1) 




1.0 


0.6420(1) 


0.63228 0.2403(1) 




1.3155 


0.5376(1) 


0.51714 0.2493(1) 




2.0 


0.3497(2) 


0.32481 0.2667(2) 





Table VII: Parameterization of the finite-size scaling function for the susceptibility Q^i^, hL/^^^^) = 
d{z) + c~^{z){hL^^/^Y^^\ determined from the results from L = 1 — 10 fm and from L = 10 — 100 

fm, for small values of the scaling variable. The last column indicates the value below which the 
parameterization is approximately valid. 



d{z) c^{z) t{z) {hLl^^l^) 



-0.5 0.30656(171) 


-1.983(888) 2.557(311) 


0.70 


-0.2 0.29607(55) 


-0.830(67) 


1.893(58) 


0.51 


-0.1 0.29311(20) 


-0.655(16) 


1.702(18) 


0.54 


0.0 0.29049(19) 


-0.535(10) 


1.532(14) 


0.56 


0.1 0.28819(46) 


-0.451(16) 


1.380(28) 


0.60 


0.2 0.28624(72) 


-0.391(18) 


1.245(36) 


0.63 


0.5 0.28207(122) 


-0.304(10) 


0.952(31) 


0.70 
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Table VIII: Parameterization of the finite-size scaling function for the susceptibility 
Q^{z, hLf^^/") = d{z) + c^{z){z, hLf^^l''Y^''\ determined from the results from L = 10 - 100 fm, for 
small values of the scaling variable. 



z 


d{z) 




t{z) 




-5.0 


0.5382(392) 


-1.567(598) 


1.702(498) 


0.52 


-4.0 


0.4887(305) 


-1.463(547) 


1.795(483) 


0.52 


-3.0 


0.4376(219) 


-1.325(460) 


1.890(444) 


0.53 


-2.0 


0.3854(138) 


-1.105(313) 


1.939(362) 


0.55 


-1.0 


0.3306(53) 


-1.134(300) 


2.162(276) 


0.53 


0.0 


0.29028(117) 


-0.4455(146) 


1.3762(368) 


0.63 


1.0 


0.28269(128) 


-0.2625(27) 


0.6135(131) 


0.68 


1.3155 


0.2891(4) 


-0.2679(3) 


0.4946(25) 


0.49 


2.0 


0.3198(95) 


-0.3111(36) 


0.3246(161) 


« 0.3 



Table IX: Parameterization of the finite-size scaling function Q^(z, /iL^^/"^) = Co^{z){hLf^^'''Y°° 
for asymptotically large values of the finite-size scaling variable, determined from the results from 
L = 1 — 10 fm and from L = 10 — 100 fm. From the infinite-volume limit, we expect that 
HQCao{z) = fxi^) agrees with the infinite- volume scaling function for the susceptibility, and that 
"Too = 5 ~ 1 is independent of z. For comparison, we list fx{z) and | — 1 in the table. 



Coo{z) 



1/6-1 Hoc^iz) f^{z) 



-0.5 0.03466(14) 
-0.2 0.03786(13) 
-0.1 0.03889(12) 
0.0 0.04230(15) 
0.1 0.04369(14) 
0.2 0.04498(15) 
0.5 0.04986(13) 



-0.8124(22) 
-0.8015(18) 
-0.7953(15) 



0.2091(8) 0.1604 
0.2284(8) 0.1831 
0.2346(7) 0.1919 



-0.8168(21) -0.7989(1) 0.2552(9) 0.2014 
-0.8142(19) ••• 0.2635(8) 0.2116 
-0.8089(19) • • • 0.2713(9) 0.2225 
-0.8056(15) • • • 0.3008(8) 0.2589 
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Table X: Parameterization for the finite-size scaling function of the susceptibility Q^{z, hL^^^^'^) = 
Coo{.z){hLf^^/'^Y°° , determined from the results from L = 10 — 100 fm, for asymptotically large values 
of the finite-size scaling variable. Prom the infinite-volume limit, we expect that Hq Coo(z) = fx{z), 
and that Tqo = 5 ~ !> independent of z. We list fx{z) and | — 1 for comparison. 



z 




'''00 


1/(5-1 HoCooiz) 




-5.0 


0.01659(14) 


-0.9616(38) 


0.1001(9) 


0.06010 


-4.0 


0.01734(12) 


-0.9307(34) 


0.1046(7) 


0.06855 


-3.0 


0.01863(11) 


-0.8950(29) 


0.1124(7) 


0.08060 


-2.0 


0.02165(12) 


-0.8656(28) 


0.1306(7) 


0.09942 


-1.0 


0.02670(10) 


-0.8265(19) 


-0.7989(1) 0.1611(6) 


0.13271 


0.0 


0.03698(7) 


-0.7933(10) 


0.2231(4) 


0.20138 


1.0 


0.05545(8) 


-0.8027(9) 


0.3345(5) 


0.31666 


1.3155 


0.05746(4) 


-0.7872(5) 


0.3466(3) 


0.33006 


2.0 


0.05137(1) 


-0.7388(1) 


0.3099(1) 


0.32481 
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